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IV 


1.  Description  of  Technology 


A  shaped  charge  jet  is  formed  by  the  detonation  of  a  high  explosive  that  surrounds  the  outside  of 
a  hollow  right  circular  cone  of  copper.  As  the  detonation  front  in  the  high  explosive  travels  from 
the  vertex  to  the  base  of  the  cone,  the  copper  is  propelled  toward  the  central  axis  of  the  cone. 

The  flow  of  copper  toward  the  central  axis  results  in  the  formation  of  a  thin  jet  of  copper  (shaped 
charge  jet)  traveling  at  a  high  velocity  along  the  central  axis  and  directed  outward  through  the 
base  of  the  cone.  This  non-molten  shaped  charge  jet  is  capable  of  penetrating  thick  annor. 

It  was  proposed  ( 1 )  that  a  large  electric  current  flowing  along  the  shaped  charge  jet  could 
vaporize  the  jet  and  make  it  harmless.  In  this  proposal,  the  electric  current  was  applied  to  the  jet 
by  two  thin  metal  plates  that  are  parallel  with  each  other  and  separated  by  some  distance.  A 
conductive  path  from  one  metal  plate  to  the  other  is  formed  after  the  jet  penetrates  each  metal 
plate.  A  current  source  that  is  connected  to  the  plates  can  now  apply  a  current  along  the  jet.  At 
the  time  of  the  proposal,  however,  a  magnetic  flux  compressor  was  the  only  type  of  current 
source  that  could  generate  a  large  enough  current  in  the  time  needed  to  vaporize  the  jet.  Since 
then,  it  has  been  shown  that  an  electric  current  can  enhance  the  inherent  hydrodynamic 
instabilities  in  the  jet  and  cause  the  jet  to  break  into  a  string  of  particles  that  cannot  penetrate  as 
much  armor  as  a  continuous  jet  (2).  It  has  also  been  shown  that  a  jet  heated  (by  an  electric 
current)  to  a  point  where  it  melts  will  break  and  expand  into  a  series  of  rings  that  have  little 
penetration  into  armor  (3).  The  current  needed  to  drive  the  hydrodynamic  instabilities  or  to  melt 
the  jet  is  far  less  than  the  current  needed  to  vaporize  the  jet.  Because  less  current  is  needed  to 
reduce  the  penetration  of  the  jet  into  armor,  it  is  possible  to  use  a  capacitor  bank  as  a  current 
source.  Indeed,  capacitor  banks  have  become  the  current  source  for  many  experiments  in 
electro-magnetic  annor  (EMA)  ( 4 ). 


2.  Experimental  Arrangement 


The  experimental  arrangement  for  many  EMA  experiments  is  illustrated  in  figure  1 .  The  parallel 
plates  are  connected  to  a  remotely  located  capacitor  bank  C  by  the  way  of  coaxial  cables,  strip 
line,  or  a  transmission  line  of  some  design.  Also  shown  in  figure  1  is  an  illustration  of  a  shaped 
charge  jet  traveling  downward  and  positioned  just  when  the  tip  of  jet,  point  A,  contacts  the 
bottom  plate  and  completes  a  conducting  path.  At  this  time,  the  current  starts  to  flow  from  the 
capacitor  bank  and  through  the  jet.  For  a  time  after  the  initial  contact,  the  current  increases  at  a 
steady  rate,  starting  from  zero  current.  Because  the  tip  of  the  jet  perforates  the  bottom  plate 
almost  immediately  after  the  current  starts,  the  current  has  little  time  to  increase  to  an  effective 
level,  and  there  is  little  time  to  drive  the  hydrodynamic  instabilities  near  the  tip  or  to  melt  the  tip. 
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In  contrast,  EMA  affects  the  portion  of  the  jet  farther  back  from  the  tip  such  as  the  portion  near 
point  B  in  figure  1 .  During  the  time  that  point  B  travels  to  the  bottom  plate,  there  is  time  for  the 
current  to  increase  to  an  effective  level,  and  there  is  time  to  drive  the  hydrodynamic  instabilities 
or  melt  the  jet.  Thus  the  rate  at  which  the  current  increases  after  the  initial  contact  must  be  as 
large  as  possible  so  that  the  length  of  the  jet  after  the  tip  that  is  not  affected  by  the  current  is 
reduced. 


The  initial  rate  of  the  current  just  after  the  initial  contact  dlo/dt  for  the  circuit  in  figure  1  is  given 
by  dl0/dt  =  V/Lt,  in  which  V  is  the  initial  voltage  on  the  capacitor  bank  and  L,  is  the  total 
inductance  of  the  circuit.  A  large  dlo/dt  would  result  by  a  large  V  and  a  small  L,.  In  this  circuit, 
the  total  inductance  Lt  includes  the  inductance  in  the  capacitors,  the  connections  to  the  plates,  the 
plates,  and  the  shaped  charge  jet.  As  a  typical  example  ( 4 ),  if  the  initial  voltage  on  the  capacitor 
bank  is  10  kV  and  if  the  dlo/dt  is  to  be  20  kA/ps,  L,  would  then  be  0.5  pH.  It  is  possible  to 
design  the  capacitors  and  the  connections  to  the  plates  so  that  their  inductances  are  a  small 
contribution  to  the  total  inductance.  The  inductance  of  the  parallel  plates  could  be  made  small  if 
their  separation  distance  were  small.  A  small  separation  distance,  however,  reduces  the 
effectiveness  of  the  armor  because  the  time  that  a  portion  of  the  jet  is  between  the  plates  and 
being  acted  upon  by  the  current  is  small.  This  reduction  in  the  effectiveness  of  the  armor  can  be 
compensated  by  a  larger  current.  In  contrast,  if  the  separation  distance  is  increased,  the  times 
that  portions  of  the  jet  are  between  the  plates  and  being  acted  upon  by  the  current  increase,  which 
results  in  less  current  to  destroy  the  jet.  Larger  separation  distances,  however,  would  increase 
the  total  inductance  of  the  circuit  resulting  in  an  increase  of  the  time  for  the  current  buildup  to 
become  effective.  During  this  buildup  time  of  the  current,  a  portion  of  the  jet  near  the  tip  is 
passing  intact  through  the  bottom  plate  in  figure  1 .  This  buildup  time  for  the  current  can  be 
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decreased  if  a  larger  charge  voltage  is  used  on  the  capacitor  bank.  Thus,  the  plate  design 
depends  on  the  power  supply  and  the  characteristics  of  the  shaped  charge  jet. 

To  design  an  EMA  system,  the  first  step  is  to  identify  a  shaped  charge  jet  to  be  defeated  and  then 
calculate  the  current  and  the  time  required  to  destroy  the  various  sections  of  the  jet.  The  next 
step  is  to  design  a  pulse  power  system  that  will  deliver  the  current  in  the  required  time,  given  the 
requirements  for  the  specific  application  of  the  EMA.  The  requirements  to  be  addressed  here  are 
the  protection  area  and  the  thickness  of  the  EMA.  Suppose  that  the  area  and  the  thickness  of 
EMA  are  given  for  some  application.  Is  it  better  to  cover  the  area  with  just  one  set  of  parallel 
plates,  or  is  it  better  to  divide  the  area  into  a  number  of  parallel  plates,  each  being  connected  to 
the  capacitor  hank  by  a  low  inductance  transmission  line?  Considering  the  inductance  of  the 
parallel  plates  in  these  two  cases  and  the  other  requirements  of  the  application  may  answer  the 
question.  Therefore,  a  method  was  developed  and  presented  here  to  calculate  the  inductance. 
Because  this  method  is  unfortunately  rather  long  and  difficult,  the  results  of  the  calculations  are 
presented  in  tables  1  through  3  for  various  plate  parameters  that  are  likely  to  be  found  in  EMA. 
With  these  tables,  which  are  presented  in  section  4,  it  is  possible  to  quickly  estimate  the 
inductance  of  a  wide  variety  of  parallel  plates  and  then  change  the  design  for  an  EMA.  Once  the 
design  of  the  EMA  has  been  narrowed,  the  inductance  of  the  plate  could  be  calculated  to  further 
refine  the  design. 


3.  Current  Distribution 


Consider  a  thin  rectangular  plate  that  has  a  current  distribution  as  shown  in  figure  2.  The  lines 
radiating  from  the  point  on  the  edge  represent  a  point  current  source  on  the  edge.  The  place 
where  the  lines  converge  to  a  point  represents  a  current  “sink”.  Each  line  is  a  flow  line  where  the 
current  density  vector  is  tangential  to  the  line.  In  EMA,  the  source  on  the  edge  of  the  plate 
would  be  where  the  plate  is  connected  to  the  capacitor  bank,  and  the  sink  would  be  where  the 
shaped  charge  jet  has  penetrated  the  plate.  The  current  distribution  on  the  opposite  parallel  plate 
is  similar:  a  current  sink  on  the  edge  of  the  plate  where  it  is  connected  to  the  return  to  the 
capacitor  bank,  and  a  current  source  where  the  shaped  charge  jet  has  penetrated  it.  The  current 
distributions  in  both  of  the  plates  are  a  source  of  a  magnetic  field  that  required  some  energy  to 
create.  This  energy  Wm  is  related  to  the  inductance  of  the  plates  by  Wm  =  LpI~/2,  in  which  Lp  is 
the  inductance  of  the  plates  and  /  is  the  current  in  the  circuit. 
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Since  different  current  distributions  in  the  plates  will  produce  different  magnetic  fields  with 
different  energies,  the  inductance  of  the  plates  will  differ.  This  means  that  the  current 
distribution  of  the  plates  must  be  known  before  the  inductance  of  the  plates  can  be  calculated. 
Unfortunately,  the  current  distributions  in  the  plates  can  be  calculated  only  after  some 
assumptions  are  made.  Each  assumption,  however,  results  in  a  different  current  distribution  and 
a  different  plate  inductance.  Furthermore,  none  of  these  assumptions  may  seem  at  first  to  be 
applicable  to  EMA.  Therefore,  these  various  assumptions  are  presented  here  and  discussed 
before  we  continue  to  the  calculation  of  the  plate  inductance. 

This  discussion  serves  as  a  tutorial  for  how  to  calculate  the  energy  stored  in  a  magnetic  field. 

The  consequences  of  using  these  assumptions  are  illustrated  by  a  common  problem  that  can  be 
solved  analytically  for  most  cases.  To  start,  consider  a  problem  as  shown  in  figure  3  that  shows 
a  right  cylinder  with  some  arbitrary  cross  section  oriented  to  be  parallel  with  the  z-axis.  The 
ends  of  the  right  cylinder  lie  on  the  z  =  1/2  plane  and  the  z  =  -1/2  plane  in  which  /  is  the  total 
length  of  the  right  cylinder.  For  all  points  in  the  cross  section,  there  is  a  current  density  vector 
J:(x,y),  parallel  to  the  z-axis  and  dependent  only  on  its  position  on  the  x-y  plane.  Now  apply  the 
general  integral  for  the  energy  stored  in  the  magnetic  field  which  is 

Wm=^\dxi  ,  (1) 
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where  the  magnetic  vector  potential  A(x)  is  produced  by  the  same  current  density  distribution 
according  to 


M, 


i 


Ax') 


A(x)  =  I  (dx'f  . 

4^l  Ax  -  x')2  +  (y  -  y')2  +  (z  -  z')2 


(2) 


This  integral  is  taken  over  the  same  entire  volume  as  equation  1 . 


Figure  3.  Right  cylinder  with  an  arbitrary  cross  section  carrying  a  current. 


When  the  length  of  the  right  cylinder  /  is  much  larger  than  the  maximum  length  across  the  cross 
section  and  when  the  current  density  is  independent  of  the  z-coordinate,  equation  1  can  be 
integrated  over  z  and  z'  to  become 


W„,  = 


Mj_ 
4  n 


\ 


dx  dy  dx'  dy'  J z  (x,  y)  J z  (V,  y’)  In 


21 


t](x-x')2  +{y-  y')2 


With  the  properties  of  the  logarithm,  this  integral  may  also  be  written  as 


W  = 


Mo  /  ln(2/)/2  n„l 


4  7t 


871 


J*  dx  dy  dx  dy'  Jz  (x,  y)  J,  (x,  y')  ln((x  -  x')2  +  (y  -  y'f ) 


in  which  /  is  the  total  current  in  the  conductor.  Now  apply  equation  4  to  a  flat  ribbon  with  a 
width  w  that  is  located  on  the  x-axis  as  shown  in  figure  4  which  gives 


W„,  = 


Mo  l  ln(2/)/2  Mol 


4  n 


4  71 


I 


dx  dx  J .  (x)  (x')  ln(jx  -  x'|)  . 


(3) 


(4) 


(5) 
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Now  an  assumption  must  be  made  that  will  result  in  a  current  distribution  Jz(x)  across  the  width 
of  the  ribbon. 


Figure  4.  Magnetic  field  lines  of  a  flat  ribbon  carrying  a  uniform  current. 

The  simplest  assumption  is  to  have  a  unifonn  current  density  across  the  ribbon  J-Jx)  =  I/w,  which 
will  result  when  the  current  is  constant  with  time  and  when  the  ribbon  is  made  of  a  material  with 
a  resistance.  After  the  double  integral  is  perfonned,  the  energy  stored  in  the  magnetic  field  for  a 
direct  current  Wf' c'  is 

m 


wd,.  =  to  1  ln(2/)  1 


4  n 


Mo  1 1" 
4  n 


(ln(w)-1.5)  . 


(6) 


The  next  assumption  is  to  have  the  ribbon  made  from  a  material  that  has  no  resistance  or  a  super 
conductor  which  changes  the  boundary  conditions  of  the  magnetic  field  on  the  ribbon.  In  this 
case,  the  magnetic  field  inside  a  super  conductor  must  be  zero,  which  is  different  from  the 
boundary  conditions  for  an  ordinary  conductor.  On  the  surface  of  an  ordinary  conductor,  the 
nonnal  component  of  the  magnetic  induction  field  B  must  be  continuous  when  we  are  going 
across  the  surface.  This  condition  also  holds  for  a  super  conductor  because  it  is  based  on  the  fact 
that  there  are  no  magnetic  monopoles.  Since  B  must  be  zero  inside  the  super  conductor,  the 
nonnal  component  of  B  on  the  surface  of  the  super  conductor  must  also  be  zero.  Therefore,  the 
magnetic  induction  field  must  be  tangential  to  the  surface  of  the  super  conductor.  For  normal 
conductors,  the  tangential  component  of  the  magnetic  field  H  must  be  continuous  when  crossing 
the  surface,  provided  that  there  are  no  surface  currents.  For  super  conductors,  however,  there  are 
surface  currents  that  can  produce  a  discontinuity  of  the  tangential  component  of  H.  The 
consequence  of  these  boundary  conditions  in  two-dimensional  geometry  is  that  the  surface  of  a 
super  conductor  is  a  field  line  of  the  magnetic  induction  field.  To  illustrate  the  differences 
between  a  normal  conductor  and  a  super  conductor,  figure  4  shows  the  magnetic  induction  field 
of  the  ribbon  with  a  unifonn  current,  and  figure  5  shows  the  magnetic  induction  field  when  the 
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ribbon  is  a  super  conductor.  For  a  uniform  current  in  figure  4,  some  of  the  field  lines  cross  the 
ribbon,  which  would  not  be  allowed  if  the  ribbon  were  a  super  conductor.  The  field  lines  in 
figure  5  were  calculated  from  a  current  distribution  that  excludes  all  the  field  lines  from  the 
ribbon.  From  other  calculations,  it  was  discovered  that  this  current  distribution  is 


which  becomes  infinite  at  the  edges  of  the  ribbon.  Substituting  this  current  distribution  into 
equation  3  gives  the  energy  stored  in  the  magnetic  field  for  a  super  conductor  W^c'  to  be 


Wl;c ■  = 


Ho  /  ln(2/)  72 


Ho  H2 
4  n 


(ln(w)  -  ln(4))  . 


Since  ln(4)  is  1.386...,  the  energy  stored  in  the  magnetic  field  for  a  ribbon  carrying  a  direct 
current  (equation  4)  is  close  to  that  of  a  super  conducting  ribbon  with  the  same  width  and  length 
and  carrying  equal  current  (equation  6).  This  is  in  spite  of  the  very  different  current  distributions 
in  the  two  ribbons. 


Finally,  it  could  be  assumed  that  there  is  an  alternating  voltage  source  connected  at  the  ends  of 
the  long  ribbon  and  operating  at  an  angular  frequency  to.  The  equation  to  be  solved  for  the 
current  distribution  across  the  ribbon  does  not  have  an  analytic  expression,  and  it  is  necessary  to 
use  numerical  methods  to  solve  the  problem.  The  results  of  the  numerical  analysis  show  that  the 
energy  stored  in  the  magnetic  field  approaches  the  value  given  by  equation  6  as  the  angular 
frequency  approaches  zero.  The  analysis  also  shows  that  the  magnetic  energy  approaches  the 
value  given  by  equation  8  as  the  angular  frequency  approaches  infinity.  This  is  the  result  of  the 
fact  that  alternating  currents  in  a  conductor  tend  to  become  more  concentrated  toward  the  surface 
as  the  angular  frequency  of  the  currents  increases.  At  very  high  frequencies,  the  currents  are 
concentrated  very  close  to  the  surface  in  the  same  manner  as  a  super  conductor.  Thus,  the 
assumption  that  the  conductor  is  a  super  conductor  gives  the  high  frequency  limit.  The 
numerical  analysis  also  shows  that  the  magnetic  energy  for  the  intermediate  angular  frequencies 
is  always  between  W^c' ,  the  lower  limit,  and  W^'c' ,  the  upper  limit. 
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Figure  5.  Magnetic  field  lines  of  a  flat  ribbon  with  a  current  distribution  given  by  equation  7. 

In  reviewing  all  of  the  previous  assumptions,  we  see  that  none  of  them  apply  directly  to  EMA 
where  the  current  is  a  transient  or  a  pulse:  it  is  not  a  direct  current  and  it  does  not  have  a  constant 
angular  frequency.  Furthermore,  the  conductors  in  EMA  are  not  super  conductors.  It  would 
seem  that  none  of  these  assumptions  could  be  applied  to  EMA  to  find  the  magnetic  energy  stored 
in  the  plates  and  then  assign  an  inductance  for  them.  The  examination  of  the  flat  ribbon, 
however,  indicates  that  the  magnetic  energy  is  not  very  sensitive  to  the  exact  current  distribution 
in  the  ribbon.  If  this  observation  holds  true  for  the  plates  in  EMA,  it  is  possible  to  find  a  good 
estimate  of  the  inductance  of  the  plates,  even  though  the  assumption  being  used  may  not  be  valid 
and  the  current  distribution  in  the  plates  may  not  be  valid.  Therefore,  it  will  now  be  assumed 
that  the  current  distribution  in  the  EMA  plates  will  be  for  direct  currents.  The  reason  for  this 
choice  is  because  the  current  distribution  for  a  direct  current  is  easily  found  with  the  use  of  any 
one  of  a  number  of  techniques.  Whatever  technique  is  used,  the  current  density  is  no  longer  an 
analytic  expression  and  the  integration  of  equations  1  and  2  must  be  done  numerically. 

Perfonning  the  integrals  in  equations  1  and  2  numerically  requires  effort.  As  an  example, 
suppose  that  the  overall  volume  is  a  cube  that  is  subdivided  into  n 3  smaller  cubes,  where  n  is  the 
number  of  smaller  cubes  along  each  edge  of  the  larger  cube.  If  the  magnetic  vector  potential  is 
to  be  estimated  at  a  corner  of  a  smaller  cube  at  (x,y,z),  the  integrand  in  equation  2  must  be 
evaluated  at  all  the  (n+ 1)  comers  (x',y',zr)  of  the  smaller  cubes  and  summed.  This  magnetic 
vector  potential  is  then  multiplied  with  the  current  density  to  become  a  value  for  the  integrand  in 
equation  5  at  this  corner.  This  procedure  must  be  repeated  for  each  corner  of  the  smaller  cubes 
at  (x,y,z).  This  means  that  the  total  number  of  times  that  the  integrand  in  equation  2  must  be 
evaluated  to  find  the  magnetic  energy  is  on  the  order  of  (n+  if .  If  the  cube  is  divided  into  100 
smaller  cubes  per  edge,  then  the  total  number  of  evaluations  of  the  integrand  is  about  10  “  which 
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is  practical  on  most  computers.  If  the  cube  is  divided  into  1000  smaller  cubes  per  edge,  then  the 
total  number  of  evaluations  of  the  integrand  is  now  about  10  ,  which  is  not  practical  even  on  a 
super  computer.  This  demonstrates  that  as  the  number  of  points  or  the  number  of  cubes  is 
increased  to  capture  the  finer  details  of  the  current  distribution  in  the  volume,  the  labor  to  find 
the  magnetic  energy  will  dramatically  increase.  Because  the  magnetic  vector  potential  and  the 
current  density  distribution  are  known  after  we  solve  a  problem,  however,  a  single  integration 
over  the  volume,  equation  1,  is  needed  to  calculate  the  energy  of  the  magnetic  field.  Finding  the 
magnetic  energy  for  two  infinitely  thin  plates,  each  with  a  current  distribution,  is  easier.  Let  the 
two  plates  be  squares,  as  an  example,  and  divide  each  plate  into  n  ,  where  n  is  the  number  of 
squares  along  each  side  of  the  plate.  The  number  of  times  that  the  integrand  in  equation  2  is  to 
be  evaluated  is  now  about  (n+ 1)4  rather  than  about  (n+l)6.  This  is  the  consequence  of  the  fact 
that  equations  1  and  2  are  now  integrals  over  the  areas  of  the  thin  plates  rather  than  integrals  over 
volumes.  Still,  a  method  was  used  here  to  estimate  the  magnetic  energy  for  a  very  large  number 
of  squares  from  the  results  using  a  more  practical  number  of  squares. 

The  result  of  using  the  procedures  based  on  finite  differences,  such  as  the  relaxation  method,  is  an 
estimate  of  the  current  densities  at  points  on  a  regular  rectangular  grid  on  the  plate.  The  plate  is 
now  divided  into  smaller  rectangular  regions  where  the  current  densities  are  given  at  each  corner. 
One  of  these  small  rectangular  regions  is  shown  in  figure  6  as  an  example.  The  origin  of  the 
coordinate  frame  for  the  rectangle  is  placed  in  the  middle  and  oriented  with  the  x-axis  parallel  to 
one  side  and  the  y-axis  parallel  to  another  side.  The  rectangle  has  a  half-length  “a”  and  a  half¬ 
height  “b”.  The  x  components  of  the  current  densities  at  the  four  corners  are  J;  ,  J'y  ,  J\'+  ,  and 

,/v  .  The  y  components  of  the  current  densities  at  the  four  corners  are  J~~ ,  Jy~  ,  Jy’+  ,  and  Jy’+  . 

The  signs  in  superscript  give  the  sign  of  the  x  and  y  directions,  respectively,  from  the  origin  to  the 
comer.  In  order  to  find  the  magnetic  vector  potential  produced  by  the  current  density  distribution 
in  this  rectangle,  the  current  density  must  be  given  at  all  points  on  the  rectangle,  but  the  current 
densities  are  known  only  at  the  corners.  Therefore,  it  is  necessary  to  assume  a  continuous  current 
distribution  that  is  specified  by  the  current  densities  at  the  comers.  Since  there  are  four  given  x 
components  of  the  current  density,  simply  assume  that 

A  O',  /)  =  « '  +  P'x  x  +  x'v'  +  S'x  x'y'  ,  (9) 

where  the  coefficients  a'x,  fi'x,  y'x,  and  S'x  are  chosen  so  that  J'x(a,b)  =  Jx'+  ,  J[(-a,b )  =  J~'+  , 
J[(a,-b)  =  J+~ ,  and  J'x(-a,-b )  =  J~~ . 


In  the  same  manner,  assume  that 

A  O',  /)  =a'v  +  /?'  x'  +  y'y'  +  S'  x'y' 


(10) 
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Figure  6.  Current  distribution  of  a  small  rectangle  on  the  plate. 


where  the  coefficients  a'y,  fi'y,  y'y,  and  dy  are  chosen  so  that  J' {a, b )  =  Jy’+ ,  J'y(-a,b )  =  J  ,+ , 
J'y(a,-b )  =  J]:  ,  and  J'y(-a,-b )  =  J ,  .  With  these  assumptions  for  the  neighboring  rectangles 


on  the  plate,  it  can  be  shown  that  the  current  density  vector  is  continuous  in  value  as  one  crosses 
the  edge  from  one  rectangle  into  a  neighboring  rectangle.  The  partial  derivatives  of  the  current 
density  vector,  however,  are  not  continuous  when  we  cross  an  edge.  In  addition,  charge  may  not 
be  conserved  within  a  rectangle,  but  it  could  be  conserved  if  Sx  =  0,  Sy  =  0,  and  f:lx  +  yy  =  0  so  that 
V  •  J  =  0  everywhere.  In  spite  of  these  limitations  of  the  assumed  current  distribution,  it  is  a 
good  approximation  of  current  distributions  that  do  conserve  the  charge  everywhere,  provided 
that  the  rectangle  is  small  enough  so  that  the  conditions  for  its  charge  conservation  are 
approximately  satisfied.  The  x-component  of  the  magnetic  vector  potential  for  this  current 
distribution  is 


u  u 

Ax  (x,  y,  z )  =  1  dx  |  dy 


,  a'x  +  j3'xx'  +  y'xy'  +  S[.x V 
tJ(x-x')2  +(y-y')2  +  z2 


and  the  y-component  of  the  magnetic  vector  potential  is 


u  u 

Ay  (x,  y,  z)  =  -^2-  J  dx  J  dy 


,  a'y  +  P'yx'  +  Y'yy'  +  S'yx'y' 
J(x  —  x')2  +  (y  —  y'Y  +  z“ 


Both  of  these  sets  of  integrals  have  the  same  form,  and  they  have  an  analytic  expression  which  is 
given  in  appendix  A.  This  rectangle  is  a  source  of  a  magnetic  field  within  which  any  current 
distribution  will  result  in  an  energy  as  stated  by  equation  1.  The  current  distribution  in  just  one 
of  the  other  rectangles  on  the  plate  will  result  in  only  a  part  of  the  total  energy  for  the  entire 
plate.  Let  the  current  distribution  in  the  other  rectangle  be 
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and 


Jx  (x,  y)  =  ax  +  fix  x  +  yxy  +  Sx  xy  , 


(13) 


Jy  (x,  y)  =ay  +  Py  x  +  Yyy  +  Sy  xy 


(14) 


The  energy  for  these  two  rectangles  is  then 

xQ  +c  ya  +d 

Mo 


H 

x0~c  ya-d 


w’v,  =  ~~  I  dx  I  dy  Ax  (x,  y,  z)(ax  +  pxx  +  yxy  +  8xxy) 


xQ+c  y0+d 


+  8?r  I  d%  I  dy  A-v  ('X’ y’  Z^a'’  +  PyX  +  yyy  +  5yXy^ 


(15) 


xo  ~c  y«  ~d 


in  which  c  is  the  half  length  and  d  is  the  half  height  of  the  other  rectangle.  It  is  understood  that 
the  two  rectangles  lie  on  parallel  planes  that  are  separated  by  the  distance  z,  and  the  center  of  the 
other  rectangle  is  located  at  (x0,  y0,  z).  It  is  also  understood  that  the  sides  of  one  rectangle  are 
parallel  or  perpendicular  to  the  sides  of  the  other.  These  integrals  are  analytic,  but  their 
expression,  also  shown  in  appendix  A,  are  long  and  complicated.  This  expression  is  applied  to 
all  possible  pairs  of  rectangles  on  the  plates  and  collected  as  a  sum  for  the  total  energy.  This 
pairing  means  that  any  two  given  rectangles  will  be  paired  together  twice:  the  first  rectangle  as  a 
source  is  paired  to  the  second  as  the  one  being  in  the  magnetic  field  of  the  first,  and  then  the 
second  rectangle  as  a  source  is  paired  to  the  first  as  the  one  being  in  the  magnetic  field  of  the 
second.  Fortunately,  the  expression  needs  to  be  evaluated  just  once  for  a  given  pair  of  rectangles 
because  it  will  give  the  same  answer  when  the  two  rectangles  are  interchanged.  Once  the 
expression  is  evaluated  for  a  pair  of  rectangles,  the  result  is  doubled  and  then  added  to  the  total. 
This  pairing  includes  each  rectangle  being  paired  with  itself,  since  each  rectangle  has  a  current 
distribution  within  a  magnetic  field  of  its  own  creation.  Here,  the  expression  is  evaluated  and 
added  to  the  total  without  doubling,  since  there  is  only  one  pairing  of  a  rectangle  with  itself. 

This  pairing  can  be  divided  into  three  groups.  The  first  group  has  pairs  of  rectangles  on  one 
plate.  The  total  energy  found  for  this  pairing  does  not  change  as  the  distance  between  the  plates 
change.  The  second  group  has  pairs  of  rectangles  on  the  opposing  plate.  The  total  energy  for 
this  second  group  is  the  same  as  the  first.  The  third  group  has  one  rectangle  on  each  separate 
plate.  The  total  energy  from  this  third  group  does  change  with  the  distance  between  the  plates. 
Once  the  total  energy,  W,  is  found,  the  inductance  of  the  plates  is  then  calculated  by  Lp  =  2.0 
W/I2 ’,  in  which  /  is  the  current  flowing  through  the  jet. 

The  explicit  expression  for  the  integral  in  equation  15  becomes  numerically  unstable  when  the 
distance  between  the  two  rectangles  becomes  large.  The  reason  for  this  instability  is  that  many 
of  the  functions  used  in  the  expression  increase  to  large  values  as  the  distance  between  the 
rectangles  increases.  The  integral,  however,  requires  these  functions  to  be  subtracted. 

Subtracting  two  very  large  numbers  that  are  nearly  equal  is  a  well-known  source  for  round-off 
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error.  At  these  large  distances,  it  is  better  and  faster  to  numerically  integrate  equation  15.  Since 
the  magnetic  vector  potential  smoothly  varies  over  the  field  rectangle  at  these  separation 
distances,  the  numerical  integration  quickly  converges.  As  the  separation  distance  is  increased 
even  more,  the  explicit  expression  for  the  magnetic  vector  also  begins  to  suffer  from  round-off 
error.  At  these  distances,  it  is  possible  to  use  a  multipole  expansion  for  the  magnetic  vector 
potential  that  is  stable  for  even  larger  distances,  but  this  was  not  needed  for  these  calculations. 

To  validate  these  calculations,  two  29-cm  by  68.5-cm  plates  were  mounted  in  parallel  by  10.0- 
cm-long  spacers  that  were  insulators.  The  plates  were  0.6  cm  thick.  One  lead  to  an  inductance 
meter  was  attached  at  the  mid  point  of  the  29-cm-long  edge  of  a  plate.  The  other  lead  to  the 
inductance  meter  was  attached  to  the  other  plate  directly  opposite  the  first  lead.  Holes  were 
drilled  through  the  plates  at  various  locations.  These  holes  allowed  for  a  brass  rod,  0.63  cm  in 
diameter,  to  be  inserted  perpendicularly  through  each  plate  and  fastened  to  make  contact  with 
each  plate.  The  inductance  was  measured  for  a  given  location  of  the  brass  rod,  and  then  the  brass 
rod  was  moved  to  another  location  for  another  measurement.  These  measurements  would  have 
included  the  inductance  of  the  plates  and  the  inductance  of  the  brass  rod,  but  the  inductance 
meter  was  adjusted  to  subtract  the  inductance  of  the  brass  rod. 

To  compare  the  calculated  inductance  with  the  measured  inductance,  we  chose  the  point  where 
the  inductance  meter  is  attached  as  the  origin  of  a  coordinate  system,  as  shown  in  figure  2,  where 
the  x-axis  is  perpendicular  to  the  edge  of  the  plate  so  that  the  x-axis  runs  along  the  line  midway 
between  the  two  68.5-cm-long  edges  of  the  plate  and  where  the  y-axis  is  parallel  to  the  29-cm- 
long  edge  of  the  plate.  Figure  7  compares  the  calculated  inductance  (the  smooth  solid  lines)  with 
the  measured  inductance  (the  squares)  for  rod  locations  along  the  x-axis  or  along  the  mid-line  of 
the  plates.  The  results  of  two  calculations  are  shown  in  figure  7.  One  calculation  assumes  that 
the  plates  are  spaced  by  10.0  cm,  but  the  plates  have  a  thickness  of  0.6  cm.  As  a  test  for  the 
effects  of  the  plate  spacing,  another  calculation  was  done  for  a  plate  spacing  of  1 1 .2  cm  which  is 
the  inside  spacing  of  the  plates  plus  the  thickness  of  both  plates.  The  comparison  shows  that  the 
calculations  can  give  a  good  estimate  of  the  plate  inductances.  In  addition,  the  calculations  and 
the  experimental  measurements  are  in  similar  agreement  when  the  rod  was  placed  at  other 
locations  away  from  the  x-axis. 
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Distance  (  cm  ) 


Figure  7.  Inductance  of  two  parallel  plates  for  hit  locations  along  the  mid  line. 


4.  Plate  Inductances 


The  inductance  of  the  plates  with  different  aspect  ratios,  different  spacings,  and  various  hit 
locations  was  calculated  and  presented  in  tables  1  through  3.  The  procedure  for  these 
calculations  started  with  choosing  the  aspect  ratio  of  the  sides  of  the  plates.  Next,  the  location 
where  the  current  is  feeding  the  plate  and  the  location  of  the  sink  or  hit  were  chosen.  The  current 
distribution  was  calculated  by  a  relaxation  method  after  the  size  of  the  grid  was  chosen.  Because 
this  current  distribution  does  not  change  as  the  spacing  between  the  plates  is  changed,  this 
current  distribution  is  used  in  the  calculation  of  the  inductance  for  a  number  of  spacings  between 
the  plates.  Keeping  the  locations  of  the  current  feed  and  the  sink  fixed,  the  current  distribution 
was  found  on  a  finer  grid  and  the  inductance  was  again  calculated  for  the  same  set  of  spacings. 
This  procedure  was  repeated  for  a  number  of  finer  grids,  ft  was  observed  that  the  inductance  for 
a  given  spacing  depended  on  the  number  of  grid  points  n  along  an  edge  of  the  plate  as  L(n)  =  Lo 
+  a/n  +  b/n2 .  An  example  of  the  inductance  converging  to  a  value  is  shown  in  figure  8  where  the 
plate  is  a  1-cm  by  1-cm  square  with  a  spacing  of  0. 1  cm.  The  current  feeds  are  at  the  mid  points 
of  two  opposing  edges  and  the  sink  is  in  the  center  of  the  square.  In  this  example,  the  inductance 
is  taken  to  be  the  extrapolated  value  at  zero  for  an  inductance  of  1 .68  nH.  All  the  inductances  in 
tables  1  through  3  resulted  in  this  extrapolation  to  zero. 
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Table  3.  1.0  cm  by  2.0  cm. 
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Figure  8.  Convergence  as  a  function  of  the  inverse  number  of 
grid  points. 

Each  table  is  the  inductance  in  nanohenries  for  plates  with  various  aspect  ratios  where  one  side  is 
always  1.0  cm.  The  spacing  between  the  plates  is  in  centimeters.  The  inductance  scales  directly 
with  the  scaling  of  the  plates.  Thus,  the  inductance  of  100-cm  by  100-cm  plates  that  are 
separated  by  10  cm  (a  scaling  factor  of  100)  gives  an  inductance  of  108.6  nH  when  the  hit 
position  is  in  the  middle  (0.50  cm,  0.0  cm),  with  the  coordinates  shown  in  figure  2.  Not  included 
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in  the  tables  is  the  fact  that  the  inductance  of  any  set  of  plates  is  zero  when  their  separation 
distance  is  zero.  With  this  fact  and  the  use  of  interpolation,  it  is  possible  to  estimate  the 
inductance  for  any  plate  spacing  spanned  by  the  tables. 

Now  consider  a  design  example.  Assume  that  a  1 .0-m  by  2.0-m  area  is  to  be  protected  by  EMA 
and  the  spacing  between  the  plates  is  0.1  m.  One  option  is  to  use  two  1.0-m  by  2.0-m  plates  and 
have  the  current  feeds  at  the  mid-points  of  the  1.0-m  edge.  Thus,  table  3  is  to  be  used  with  a 
scale  factor  of  100.  The  hit  location  with  the  highest  inductance  is  always  at  the  farthest  point 
from  the  current  feeds,  x  =  2.0  cm  y  =  0.5  cm,  and  z  =  0. 1  cm  in  table  3.  Multiplying  the 
inductance  from  table  3  (3.7  nH)  by  the  scale  factor  gives  a  plate  inductance  of  370  nH.  This 
inductance  can  be  reduced  by  the  use  of  two  sets  of  parallel  plates  that  are  1.0  m  by  1.0  m,  each 
with  its  own  current  feeds.  Thus,  table  2  is  now  to  be  used  with  a  scale  factor  of  100.  Again, 
using  the  farthest  point,  x  =  1.0  cm,  y  =  0.5  cm,  and  z  =  0.1  cm,  for  the  maximum  inductance 
(2.5  nH),  the  maximum  plate  inductance  is  reduced  to  250  nH.  This  reduction  of  the  inductance 
must  be  weighed  against  the  fact  that  another  current  feed  must  be  added  to  the  system.  A  final 
option  is  to  again  use  1.0-m  by  2.0-m  plates  but  have  the  current  feed  at  the  mid-point  of  the 
2.0-m  edge.  This  gives  the  aspect  ratio  of  0.5  for  the  plates  (table  1)  and  a  scale  factor  of  200. 
The  location  of  the  maximum  inductance  would  then  be  when  x  =  0.5  cm  and  y  =  0.5  cm  but 
with  z  =  0.05  cm  for  the  spacing,  which  is  not  tabulated.  Using  the  fact  that  the  plate  inductance 
is  zero  when  z  =  0.0  cm  and  observing  that  the  inductance  increases  roughly  as  the  square  root  of 
z  when  z  is  near  zero,  the  inductance  for  z  =  0.05  cm  is  approximately  0.707  times  2.25 1  nH,  the 
inductance  for  z  =  0. 1  cm,  for  an  estimated  value  of  1 .59  nH.  The  calculated  inductance  for  this 
spacing  is  1 .47  nH,  which  shows  that  this  interpolation  gives  a  reasonable  estimate.  This  gives 
the  estimated  maximum  inductance  for  this  option  as  3 18  nH,  while  the  calculated  maximum 
inductance  is  294  nH. 


5.  Conclusion 


A  procedure  to  estimate  the  inductance  of  two  parallel  plates  that  have  a  current  distribution  on 
them  has  been  presented.  Because  the  equations  used  in  this  procedure  are  long  and 
complicated,  tables  of  inductances  for  a  set  of  parallel  plates  are  presented.  By  interpolating  the 
inductances  in  these  tables  and  using  a  scaling  law,  it  is  possible  to  estimate  the  inductance  of 
any  set  of  parallel  plates  to  within  about  10%,  which  is  a  sufficient  tolerance  for  most  cases. 
These  equations  could  be  extended  to  find  the  inductance  of  more  complicated  structures,  but 
they  may  become  even  more  complicated. 
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Appendix  A.  Integrals  for  the  Inductance  Calculation 


The  integration  of  equation  1  for  the  energy  stored  in  the  magnetic  field  starts  when  the 
expression  for  the  magnetic  vector  potential  in  equations  1 1  and  12  is  found.  Since  equations  1 1 
and  12  become  identical  when  the  x  and  y  subscripts  of  the  coefficients  were  dropped  and  when 
the  limits  of  integration  are  ignored,  only  equation  1 1  is  integrated  here  with  its  x  subscripts 
dropped. 

=  f  W\#  f+V+ry+W  ,  .  ,1-A) 

4n  J  J  a/(x-x)  +(>>- v)  +Z- 


(1-A) 


As  a  reminder,  the  rectangle  is  on  the  x-y  plane  with  its  center  at  the  origin  of  the  coordinate 
system.  The  normal  vector  of  the  rectangle  is  in  the  positive  z  direction.  The  rectangle  is 
oriented  with  one  set  of  its  parallel  sides  perpendicular  to  the  x-axis.  The  half  length  of  the 
rectangle  along  the  x'  axis  is  a  ’,  and  the  half  width  of  the  rectangle  along  the  y  axis  is  b  ’.  The 
integration  starts  when  the  variables  are  changed  to  x'  =  x  -  u  and  y'  =  y  -  v.  After  these  new 
variables  are  substituted  into  the  equation  and  after  the  like  powers  of  u  and  v  are  collected  in  the 
numerator  of  the  integrand,  the  equation  becomes 


A(x,y,z)  = 


ju0(a'  +  /3'x  +  y'y  +  S'xy) 


x-a  y-b' 


AtiM  L„j\ 

x+a  y+b' 
x-a'  y~b' 

-  I* du  \dv- 

J  J  , 

x+a'  y+b' 
x-a'  y~b' 

MAA  \du  L 

J  J  A 


x+a'  y+b' 
x-a'  y-b' 


u2+v2+z2 


2  .  2  .  2 

\u  +  V  +z 


/  2  .  2  .  2 

\u  +  V  +z 


-  I du  I  < 

x+a'  y+b' 


/  2  .  2  .  2 

In  +v  +z 


(2-A) 


The  four  double  integrals  can  be  written  in  a  general  fonn: 


,(k,v,z)  =  \du\‘ 


u2+v2  +  z 2 


(3 -A) 


in  which  the  subscripts  i  and  j  are  the  powers  of  u  and  v,  respectively. 


Letting  r  =  y/u2  +v2  +  z2  ,  the  expression  for  the  integrals  used  in  equation  2-A  is 


19 


and 


/00(w,v,  z)=[du[ 

Il0(u,v,z)  =  M 
I0l(u,v,z) 
Iu(u,v,z)  = 


=M 


dv 

dv 

dv 

dv 


I  2  .  2 

V  U  +  V  +  z 


a/m2  +v2  +z2 


=  u  ln(v  +  r)  +  v  In (u  +  r)-z  arctanj 
vr  +  (u2  +  z2)  ln(v  +  r) 


b  uv^ 


yzrj 


ur  +  (v  +  z  )  ln(fi  +  r ) 


a hr  +v2  +z2 
uv 

a/m2  +v2  +z2 


With  these  definitions  of  the  integrals,  equation  2- A  can  be  written  as 

Af  x  juja' +  fix  +  y'y  +  S'xy)  r 

-f’ z)  = - 7- - Ko(«  A,  z)L,,,+, 


4^ 


L  l,uv  lx+a,y+b 

4^  L  U>1V  7Jx+a,y+fe 


Mo  S' 

-^—Viuyu,v^)ix+a,y+h,  . 


[lu(u,v,z)l~+aj 


y~b' 


(4-A) 


The  brackets  mean  that  the  expression  inside  is  evaluated  at  the  limits  of  integration.  Normal 
notation  would  have  only  one  number  in  the  superscript  and  one  number  in  the  subscript,  but 
here  there  are  two.  In  this  case,  the  brackets  mean  the  following: 

\_F(u,v,z)^+a^i’^+^y  =  F(x  -  a',y  -  b',z)  -  F(x  +  a',y  -  b',z)  -  F(x  -  a',y  +  b',z)  +  F(x  +  a',y  +  b',z)  .  (5-A) 


Note  that 


\1  x-a',y-b'  T  xlx+a',y+b'  T  r-»/  '~a',y+b'  f  n/  ~\x+a',y-b' 

[F(u,v,z)\x+a,y+b,  =  [F(u,v,z)\x_a,^b,  =-[F{u,v,z)l+a,y_b,=-[F{u,v,zl_a,y+b,. 


Because  of  round-off  error,  equation  4-A  should  be  used  with  caution  when  we  are  calculating 
the  magnetic  vector  potential  at  a  field  point  (x,y,z)  far  from  the  rectangle.  The  source  of  this 
round-off  error  comes  from  the  subtracted  terms  in  the  definition  of  the  brackets  equation  5-A 
such  as  F(x-a  \y~b  \z)  -  F(x+a  ’,y~b  \z),  which  means  that  the  function  Iij(u,  v,z),  as  an 
example,  is  to  be  evaluated  at  points  (x-a’,y-b’,z)  and  (x+a’,y-b’,z)  and  then  subtracted,  that  is, 
Iij(x-a  ’,y~b  \z)  -  Iij(x+a  ’,y~b  \z).  Unfortunately,  the  value  for  Iu(x-a  ’,y~b  \z)  is  large  and 
nearly  equal  to  Iyi(x+a  ’,y~b  \z)  at  field  points  far  from  the  rectangle.  Subtracting  two  very  large 
numbers  that  are  nearly  equal  is  a  well-known  source  of  round-off  error  in  computers. 
Fortunately  for  the  far  field  points,  the  magnetic  vector  potential  can  be  accurately  calculated 
with  the  use  of  its  multipole  expansion  which  does  not  suffer  from  round-off  error. 
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To  find  the  energy  stored  in  the  magnetic  field,  equation  4-A  is  multiplied  by  the  current 
distribution  in  the  second  plate,  J{x, y)  =  a  +  /3x  +  yy  +  8xy,  and  integrated  over  the  region  of 

the  second  plate.  The  center  of  the  second  plate  is  located  at  the  point  (x0,  y0,  z)  and  has  a  half 
length  a  and  a  half  height  b.  This  multiplication  would  result  in  an  awkward  expression.  To 
avoid  this  awkward  expression,  the  multiplication  and  integration  of  each  term  in  the  above 
function  J(x,y)  are  perfonned  separately.  Thus  let 

xQ+a  yQ  +b 

W00(xo,yo,z)  =  —  \dx  J dy  A(x,y,z)  , 

xa  -a  yQ  —b 
n  xa  +a  ya  +b 

w\,0  (X0’T0?Z)  =  —  \dx  jdyxA(x,y,z), 

xq  -a  yQ  —b 
xQ+a  ya+b 

K,1  (xo’y0’z)  =  Y  \dx  \dy  y  A(x,y,z), 

xQ-a  y0~b 


The  integration  of  W \  j/x0,y0,z)  is  given  in  detail,  which  starts  with  the  substitution  of  equation 
2-A  into  the  above  integral  for  W\fi(x0,y0,z),  but  only  the  final  results  for  the  rest  are  presented. 


xQ+a  ya+b 


Wyo  (xo  ,y0’z)  =  ^dx^dyx  (a!  +  fix + y'y  +  S'xy)[l 0  0  (u,  v,  z) J^/^, 

xn—a  yn—b 


toP 

8  n 

toP 

8  n 


xQ  +a  yQ  +b 


M  dy  x  (/?'  +  8'y)[t  0(u,  v,  z)J+“ 


y-b' 

y+b' 


xa  -a  yQ  —b 
xQ+a  yQ+b 


\dx\  dyxiy'  +  Sfyl 


0~a  y0-b 
xQ  +a  yQ  +b 


-  \dx  J^x[7u(m’v’z)K;^ 


(6-A) 


XQ  —a  yQ  —b 


All  the  double  integrals  could  be  written  in  a  short  form  as 

x0+d  ya+b 


cyn.m  I  i  «  ffl  r/  \  ]x-a'  ■ 

s‘j  =  \dx  \dyx  y  [/,  ;(u,v, z)^. 


y-b' 

y+b' 


(7-A) 


x0  -a  yQ  -b 


where  the  superscripts  n  and  m  are  the  powers  of  x  and  y,  respectively,  and  the  subscripts  are  the 
subscripts  of  the  double  integral  defined  by  equation  3-A.  With  this  notation, 
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y. .  *) = <“'  s« + P‘ C + /  s0“  +  <?'  si  i ) 

67t 

Ao  ^  t  Q1  oO’O  i  ST'  C^'*  i  *,r  oO.O  i  £r  ol>0  £r  o0’0\ 

-— — (p  >1,0  +s  >1,0  +r  >0,1  +d  Sox  -5  Sxx  ) 

8  n 


(8-A) 


wu(x<,,y„, z)  =  SJJ  +  P  s0,'°  +  /  sio  +  £'soi) 

07T 

M o  P  /  Qt  ol.O  .  or  rrl.l  .  r  rrl.O  .  or  a 2,0  or  ol,0\ 

— - — (p  >i,0  +  5  >i,o  +  r  ->0,1  +  O  Sox  -  8  Sxx ) 

8  n 


(9-A) 


KP»,y».z)=fz(a'  «s + /S'  sii + /  sg + <5'  ® 

o;r 

Ao  y  {  Q1  o0,l  I  0^.2  I  „/  o0,l  i  oM  O'  o0,l\ 

- —  IA  *>1.0  +S  >1,0  +  ^  >0.1  +  ^  >0,1  -<>  >1.1  ) 

8;r 


(10- A) 


and 


Wu(x„y„,z)  =  ^(a'  SJJ  +  ^Sg  +  r’S“+*’  Sg> 
8;r 

- ^  (/?'  Sg + S'  s;~  +  /  si: + S'  sfi  -  S'  sg  > 

8;r 


(11-A) 


Thus,  there  are  a  total  of  24  different  double  integrals  to  be  evaluated.  The  procedure  for 
evaluating  all  these  double  integrals  starts  with  the  removal  of  the  bracket  notation  in  equation 
7-A  and  the  terms  written  explicitly. 


xQ+a  ya+b 


rin,m  _ 

*uj  ~ 


=  1*1 


dy  x" ym  It  . (. x-a',y-b z ) 


xQ  -a  yQ  -b 
xQ  +a  ya  +b 


xQ  —a  yQ  -b 
xQ+a  yQ+b 

xa  —a  yQ  -b 
xQ+a  y0+b 

+  JAjsvxV"/,>  +  a',v  +  S,z). 

xa -a  yQ  -b 


The  next  step  is  to  change  the  variables  in  each  term: 

x0+a-a'  y0+b-b' 


rin,m  _ 

Aj  - 


=  J dp  J ’ dq  (p  +  a'y(q  +  b')m  ILj(p,q,z) 

xQ-a-a'  y0—b-b' 
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xQ+a+a'  y0+b-b' 

[dp  J ’  dq(p-  a')"  ( q  +  b')m  /,  f{p,q,z) 


xQ-a+a'  y0-b—b' 
xQ+a-a'  ya+b+b' 


\dp  j 


dq  ( p  +  a')’\q-b')m  b  Ap,q,z ) 


x0—a-a'  y0-b+b' 
xQ+a+a'  ya+b+b' 


fdp  J 


dq  (p ~ a')" (q  —  b')m  ILj(p,q,z)  . 


(12- A) 


x0-a+a'  y0-b+b' 

The  term  in  front  of  each  Ijj(p,q,z)  can  be  multiplied  into  a  sum  of  tenns  such  as 

k=n,l=m 

(P  +  a')n(q  +  b')m  =  ^  '  Ck  lpkq' .  Thus  the  double  integrals  in  equation  12-A  now  have  the  form 


i-=0,/=0 


H?j  = 


-=I4 


dq  pkq'lij(p,q,z) 


(13-A) 


whose  expressions  can  now  be  determined.  With  the  notation  defined  by  equation  13-A,  the 
integral  Sff  is  given  as  an  example.  Thus, 

xQ+a-a'  y0+b-b' 

S?j  =  J  dp  ^dq  ( p 2  +  lap  +  (a')2)  ILj(p,q,z ) 

x0-a-a'  y„  -b-b' 
x0+a+a'  yQ+b-b' 


[dp  J ’ dq{p2  -  lap  +  (a')2)  ItJ(p,  q,  z) 


xQ—a+a'  y0—b—b' 
x0+a-a'  y0+b+b' 


[dp  [dq  {pz  +  lap  +  (a')z)  IUJ(p, q, z) 


x0-a-a'  yQ—  b+b' 
x0+a+a'  yQ+b+b' 


\dp  J‘ 


dq(p2  -lap  +  (a  )2)  b  ( p,q,z ). 


IJ 


xQ-a+a'  y0-b+b' 


After  the  definition  for  bracketed  terms  and  the  definition  for  the  double  integrals  are  used, 
equation  1 1-A,  this  expression  becomes 


s?:‘ .  [»■« + 2a'H]j + -  k?  -  t 

-  k; + uh»  *  + k?  -  + <» 


+a+a',^0  +b—b' 
j  4 xn-a+a' ,ya-b-b' 


+a+a',y0  +b+b' 
:0  -a+a',y0  -b+b' 


(14- A) 
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where  the  functions  H"’”\p,q,z)  were  abbreviated  to  .  Note  that  equation  12-A  holds  for 
any  values  for  i  and  j.  Thus  there  are  a  total  of  nine  unique  forms  for  S"'m  in  equations  6- A 
through  9- A.  The  other  eight  forms  are 


<-,0,0  _  [zj-O.oK  +a-a',ya+b-b 
~  L ni,j  \Xo-a-a',y0-b-b: 


Q  +a+a'  ,ya  +b-b 
-a+a'  ,yQ  —b—b' 


4x0  -a-a',y0  -  b—b'  4x0  -a+a',yQ  -  b—b'  L  z> 

y;  -  k; + -  k”  - 


’  -  \]j0'0'\°+a-a'’y°+b+b'  +  [//0-0]V° 
L  '  ./  ^-a-a'^-A+y  L  ‘J  J.V0 


k 


k=k;‘+*A;“];; 


,  +<z+a',}>0  +b—b' 

-uo  ~  -  -o  -ci+a' ,y0-b-b' 

a'H0'0}o+a-a’y°+b+b  +  Ih1’0  +  a'H°’°]‘°+a+a'’y°+b+b 

l’J  *x0  -a-a',y0  -b+b'  L  l’J  l’J  *x0  -a+a',y0  -b+b' 

+a-a',y0  +b—b' 

-a-a'  ,yQ  —b—b' 


+a+a',y0  +b+b' 
a+a',yQ  -b+b' 


b'  _  \h*1  +  ij'fjO,olFo+a+a'’yo+b-b' 
-  -Xo-u-u  ,Xo-u-b'  1  ’’j  *’j  -1 lc0-a+a',y0-b-b' 

-  \h*)  -  b'Hf’f  1*"  +a~a'’yo  +b+b'  +  \H) 

L  l’J  l’J  -lx0  -a-a  ,yQ  -b+b  L  1 


70,1 

i,j 


a-a'  ,yQ  —b—b' 
+a+a'  ,ya  +b-b' 
-a+a',yQ  —b—b' 
+  a  -  a' ,y 0  +b  +  b' 
a-a' ,ya  -b  +  b' 


+ 

A!=k 

-ky- 


+  b'H*’*  ]v°  +a~a’y° +b+b' 

’’j  -tx0—a—a,y0—D+o  ■-  l’J  *x0-a+a',y0-b+b' 

s;i = k; 

-  [H"+a'H?j -b'H^-a'b'Hf&+a-a;y^l 

L  1>J  l’J  l’J  l’J  Jx0  -a-a  ,yQ  -b+b 

+  \h))  -a'//,0'* -b'H)'* +d  b'H*'* f+a-a'^+b+b) 

L  l’J  l’J  l’J  l’J  -*x0 -a+a ,yQ-b+b 

S2)  =  \ll2)  +  2  a' H]’)  +  (, a'fH *)  +  b'H2’*  +  2  d  b'H)*  +  (df  H*’*^+a~a';y°+b~b', 

l’J  L  l’J  lH  ^  7  l’J  l’J  l’J  ^  7  1>J  *x0-a-a  ,ya-b-b 

-  [ll2j  -  2a'H))  +  (a')2#  +  b'tf2;°  -  2d  b'H);*  +  (dfH*'*  \ 

-  \h )j  +  2a’H))  +  (a')2H°j  -b'Hff  -  2a' b'H ))  +  ( ,/ ' )  // ^' 1 

k>’  -  2»ws+(«v«s  22 

” 0  2b’ H*2  +  (b')2H°-°Y+a-a'’y°+b-b'  -  \h°?  +  2b'H*1 

1>J  l’J  4x0  -a-a  ,yQ  -b-b  L  lH  *»' 

2 b'H0'1  +  (b')2H°’ *\°+a-a''y°+b+b'  +  [//0.2  _ 

v  7  l»J  *x0-a-a  ,yQ-b+b  L  N7  _.o 

S''2  =  \H)'2  +  2  b'H)')  +  {b')2  H)’°.  +a'Hj’f  +  2a' b'H*')  +  {b'f  H^+a-a'^+b-b\ 

lJ  L  l’J  l’J  v  ’  '•!  l’J  1>J  v  ’  l’J  ix0-a-a  ,ya-b-b 

-  [ H)'2  +  2b' H))  +  (b'  fll)) - a'H ,°f  -  2 a' b'H*)  +  (b'f  H*’*]°+a+a';y°+b~b' 

L  1>J  hj  v  7  hj  i,J  hj  v  7  hj  ■iXo-a+a  ,yQ-b-b 

-\h)'2  -2b' H))  +  (b'f  H)'*  +  a'H*'2  -2a' b'H*)  +  (b'f  H*’*]°+a^;y°+b+b', 

L  hj  i,J  ^  '  i,J  i,J  hj  ^  S  i,j  J (xo-a-a 

+  [//'f  -  2 ////'  '  +  (/f)2//1'1  -a'H*'2  +  2a' b'H*)  +  (h')2^,0’“T°' +a+a';y°+b+b' 

L  x  y  l’J  lH  l’J  x  7  l’J  J xQ-a+a  ,yQ-b+b 

~H2;2  +  2  a'H)’2  +  (a1)2  H*;2  +  2b' II2)  +  4  a'  b'H)’)  + 

2  (a')2  b'H*)  +  ( b'fH :;*  +  2a' (b'f  H’*  +  {a’  b')2  H*;* 

~H2’2  -  2 a'H)’2  +  (a1)2 H*j  +  2b' II2)  -  4 a  b'H))  + 

2  (a'fb'Hfj  +  (b'fH2'*  -  2d  (b'f  H)’*  +  {d  b')2  H*;*  y 


u  -a+a' ,yQ-b+b' 

,  fU'\2  U°,o}C°+a+a''y°+b-b' 
j+{b)  H,  j  lo_a+a,  yo_b_b, 

2  bH°)  +  (a'fH°:°l\ 


a  +  a'  ,yQ  —b—b' 
c0  +a+a'  ,y0  +b+b' 
lxn  -a  +  a',y0  —b+b' 


k = 


-a-a' ,yQ  -b-b' 
\xa  +a+a',yQ  +b—b' 
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II;;:  +  2a’Il';j  +  (, a')2  H °;2  -2b'H;j  -  Ad  b'H )f.  - 
2(a’yb'H(;j  +  (b'fHff  +  2a\b'fH)fiJ  +  (, a'b'fH ®-° 


J  . 

H2;2  - 2a'H'j  +  (a)2 H°j  -2b' H2;)  +  Ad b'H)\  - 
2 (dfb'H*;)  +  ( b')2H 2;°  -  2d(b'fH)'°j  +  (d  b'fH °;° 


x0  +a-a',yQ  +b+b' 


x0  —a— a' ,yQ  -b+b' 
x0  +a+a',y0  +b+b' 


x0  ~a+a\y0  ~b+b' 


These  forms  or  expressions  for  S"’™  are  substituted  into  equations  6- A  through  9- A  with  the 


appropriate  values  for  the  indexes.  After  this  substitution,  there  are  a  number  of  Hd  that  appear 


more  than  once.  After  we  take  into  account  the  duplicate  integrals,  there  are  25  unique  Hf'j , 


which  are  given  below  where  r  =  pJ  +  q2  +  z2  . 


H°’°  =  — (2 z2  -  p1  -q2)  +  —  (p2  -z2)ln(q  +  r)  +  —  (q2  -  z2)\n(p  +  r)  -  pqz  arctanl 


V  rz  J 


H 0.0  =  dd  z 2  ~cl2  -  6 P2)  +  ^d~ ln(^  + 

4o  j 

+  -l(6<,y-3Z4-12pV+,4  +  12p¥)ln(P  +  r)-f(r+3/)arc,an| 


f  p£ 
yrz  ) 


H02  =  Z: (3^2  - p2  ~ 6 q2)  +  —  (6/7  V  - 3 z4  - 12 q2z2  +  p4  + 12 p2q2)\n(q  +  r) 


1 


48 


48 


Pd 


zp ,  2 


+  In (p  +  r) - (z  +  3</  )arctan| 

3  6 


( p£ 
V  rz  J 


Hoi  =^p(z2  - P2  -q2)  +  ^-(5z2  +  p2  +\0q2) In (q  +  r)  +  (5z2  + 1  Op2  +  q2 ) In (p  +  r) 

60  60  60 


- (z  +5 p  z~  +5q"z~  +  15/rg~)arctan| 

60 


(  pcf 
V  rz  J 


H 


0,0  _  Kj_  __2  o  „  2  c  _  2  \  ,  1  /  /f  „  2  _  2  ,  4 

0,1 

pq 


24  (3/7-2g-5z)+  —  (6p  z  +  p  -  3z  )ln(  q  +  r) 


(3  z 2  +  g  2 )  ln(  /?  +  r  )  - 


^z 


-  arctan 


f— 1 

V  >'Z  J 


*i?=^(2V  -3r/2-7z2)  +  -^-(5z2  +/72)ln(^  +  r) 


9 


+  240  +  60/?  z  +  3</  +  20/?“g  )ln(/?  +  r)-  —  (z  +  5/?“)arctan 


f— 1 

V  rz  J 


=  TTA (9 pV  -  8z4  - 1 6</2  z2  -  8</4  +  9p2q2  +  2p^)  +  I-(q2  +  z2y  ln(p  +  r) 


,2  _2  o_4 


22 


2  „2 


,4\  ,  P  ,  „2x2- 


120 


8 
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Hi\  =^r(V  +i6pV -3q‘  +  16+z!  -6,V  -3z4) 


288 


1 


+  I^(2z  +  6q  z  +\2p  z  +  6q  z  +Up  q  z  +2q  +  l2p-q)Hp+r) 


=  -^(3C  -  2p2-5z2)  +  £3-(p2  +  3zJ)ln (q  +  r) 
24  6 

+  -^j-(g4  +  6g2z2  -3z4)ln(/;  +  7")-^-arctan| 


f  p£ 

V  rz  J 


H w  =  ^  (22g2  -  3/  -  7z2 )  +  (1 5z4  +  60 q2z2  + 1 0p2z2  +  20 p2q2  +  3p4)  In (q  +  r ) 


240 


240 


+  30  +  5z^ln(p  +  r^3Q<<q2+  5z^  arctanl 


f  p£ 
V  rz  J 


HI'q  =  j^-(9g2z2  -^^p^z2  +  2  q4  +9  p2q2  -8  p4)  +  -^(z2  +u2)2  ln(^  +  r) 


H\’l  =  (4 q4  +  \6p2q2  -3 p4  -3z4  +  \6q2z2  -6 p2z2) 


1 


288 


H - (z  +  6 q"z  +3 p  z  +  \2p-q  z  +3 p  z  +6 p  q  +  p  )ln(q  +  r ) 

96 


Kf  =  (1 8z2  +  lq2  +  Ip2 )  +  (1 5z4  + 1 0p2z2  +  3p4)  In (q  +  r) 


120 


120 


+ - (15z4  +  10g2z2  +3qi)h\(p  +  r ) - arctan 

120  15  1 


(  pq4 

V  rz  J 


//,CV  =  (33z4  +  66  q2z2  +  26  p2z2  +33q4  +  26p2q2  +Sp4)  +  —(z2  +  q2f  In  (p  +  r) 


720 


48 


r+0 _  V  /z)z)_4  ,  o/r_2_2  ,  /r/r„2_2  ,  o  _4  ,  '->zr__2_2  ,  'n_4\  ,  ^  /  _2 


// ’,u  =  -222—  (33z4  +  26 q2 zl  +  66 p' zl  +  Sq4  +26 p2q2  +33 /?4)  +  —  (z2  +  /7°)3  ln(<7  +  7') 


720 


48 


HU  = 


105 


H™  =  —  (9q2z2  -  8z4  +  4/rz2  +  2<?4  -  p2*2  - 18/)  +  £ (z4  +  3p4) \n{q  +  r) 


180 

+  ^(r-z2)ln(^  +  7-)-^arctanM 
6  3  \  zr ) 


12 


26 


772;,0  =  (99 z4  +  78 q2z2  -  27 p2z 2  +  24^4  - \2p 2q2  +  234 /) 

+  — (z6  +  9/?4z2  +  2/?6)ln(g  +  r)  +  ^— ^-(3z2  +  c/2)ln(/;  +  r)  -  ^  -  arctan  — 

72  18  9  vzr7 

77 2’o  =  ^  (i  6^2z2  -  3z4  +  Vz2  +  V  -  2^  V  -  6/) 

i  3  3 

+  j^-(z6  +  6g2z4  +  9p4z2  +  18 p4q2  +  2 //’ )  ln(  v  +  r)  +  ln(  p  +  r) 

-^(z2+3r)arctanf^l 
18  v  zr  y 

77^  =  -2^-(16/rz2  -3 z4  +  3g2z2  +  4/?4  -2 p2 q2  - 6c/4 )  +  ^  9  ln(v  +  r) 

+ - (z6  +  6p2z4  +9q4z2  +  18p2q4  +  2q6)ln(u  +  r)-^—{z2  +  3/r)arctan  — 

144  18  ^  zr  J 

<02  =  (9/>  V  -  Bz4  +  4q2z2  +  2 p4  -  p2q2  - 1 8<?4)  +  ^ (^4  +  3<?4) Hp  +  r) 

+  ^(/-z2)ln(g  +  .)-^arctanf^l 
6  3  V zr  7 

<  =  -j^(33z4  +  26Plzl  ~  V*2  +  8 p4  ~  4p2q 2  +  78 q4)  +  ^|-(3z2  +  P2)Mq  +  r) 

+  — (z6  +  9g4z2  +  2q6)\n(p  +  r)-9  —  arctan  -^L 
72  9  l  zr  J 


=  —  (1 6z6  -  8 q2z4  -  8 p2z4  + 1 1 1  q4z2  +  12 p2q2z2  +  111  p  V  +  30^ 


,4  2  .  on  J 


3  3 

-15 p2q4  -15  p4q2  +30 p6)  +  —  (z4  +2p4)\n(q  +  r)  +  —  (z4  +2q4)\n(p  +  r) 

36  36 

-£77arctan('i>£| 


77-  =  (1 6z6  -  8q: 2z4  +  48/rz4  + 1 1 1  q4z2  -16p2q2z2  +  4Sp4z2 +30q6 

+ 1 1  \p2q4  -  %p4q2  + 1 6p6)  +  (z2  +  p2)2  In (q  +  r) 


27 


Hlj  =  (1 6z6  -  8p2z4  +  48q2z4  + 1 1  lp4z2  - 1 6p 2q2z2  +  48?  V  +  30 p6 

n3 

+ 1 11  p4q2  -8p2q4  +  \6qb)  +  —  {z2  +  q2)2  In (p  +  r ) 
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